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Abstract 

The Wilson discretization of the dimensionally reduced supersym- 
metric Yang-Mills theory is constructed. This gives a lattice version of 
the matrix model of M-theory. An SU(2) model is studied numerically 
in the quenched approximation for D=4. The system shows canonical 
scaling in the continuum limit. A clear signal for a prototype of the 
"black hole to strings" phase transition is found. The pseudocritical 
temperature is determined and the temperature dependence of the to- 
tal size of the system is measured in both phases. Further applications 
are outlined. 
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1. Introduction and lattice formulation. The matrix model of M-theory 
provides a qualitative description of many essential properties of the final 
unifying theory |]J. The interest in studying quantitatively its properties 
is particularly enhanced due to a rich spectrum of predictions from super- 
gravity, for example the thermodynamical properties extracted from black 
hole solutions, for the system of DO branes described by this SYM quantum 
mechanics 0, [|. 

The complete solution of the above matrix model is not known even 
though it is much simpler than a conventional field theory. In this letter we 
construct the Wilson discretization of the model, and study its yet simpler 
version with the well developed lattice methods. In particular, we find the 
onset of a prototype black hole to strings phase transition and determine 
quantitatively some of its properties. Of course the true black hole phase of 
the full DO brane system is much more complex. It is therefore quite appeal- 
ing that the simplified model considered in this exploratory study reveals an 
important part of the phase structure. 

Many other problems can be attacked within the present approach open- 
ing a new area of exciting applications. To our knowledge, this is the first 
study of the M-theory related quantum mechanics on a lattice. The zero di- 
mensional system has been considered recently in an attempt to understand 
the compactification of higher dimensions |5|. Monte Carlo and analytical 
study of the relevant SU(N) integrals have been reported in ||. 

We begin with the Banks, Fishier, Susskind and Shenker (BFSS) proposi- 
tion to use the dimensionally reduced SUSY YM theory in D=10 dimensions 
as a model for the relevant degrees of freedom of M-theory. For general D 
the action reads |7J 

S = Jdt QTri^) 2 + * a (t)V* a (t)} . (1) 

In the process of dimensional reduction all fields are assumed to be indepen- 
dent of the space variables x i: i = 1 . . . D — 1. Consequently all space deriva- 
tives in the field tensor F^ v and in the Dirac operator T> vanish (di — 0), 
and (|l|) describes supersymmetric quantum mechanics of D — 1 bosons and 
their fermionic partners. The temporal components of the gauge fields are 
nondynamical and serve to impose Gauss law constraints. The original D 
dimensional theory is supersymmetric at the classical level only in D=2, 4, 6 
and 10 dimensions, where appropriate (Majorana, Weyl or both) conditions 
are imposed |§. Fermionic fields ^ a (t) belong to the adjoint representation 
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of the gauge group SU(N), a = 1 . . . N 2 — 1. Finally, the BFSS proposal 
requires N — > oo since this variable corresponds to the eleventh component 
of the momentum in the infinite momentum frame where the original theory 
is considered. 

We propose to study the above system with methods of the Lattice Field 
Theory. To this end consider D dimensional hypercubic lattice Ni x . . . x Nd 
reduced in all space directions to N{ = l,i = 1 . . . D— 1. Gauge and fermionic 
variables are assigned to links and sites of the new elongated lattice in the 
standard manner. The gauge part of the action reads 

N t i 

Sg = ~PY: E^(Tr^(m)), (2) 

m=l p>v 

with 

P = 2N/a 3 g 2 , (3) 

and Ufj, u (m) = Ul(m)UUm + u)U v (m + fi)U^(m), U^ijn) = exp (iagA^am)), 
where a denotes the lattice constant and g is the gauge coupling in one 
dimension. The integer time coordinate along the lattice is m. Periodic 
boundary conditions U^m + v) — U^m), v = 1 ... D — 1, guarantee that 
Wilson plaquettes tend, in the classical continuum limit, to the appro- 
priate components with space derivatives absent. In this formulation the 
projection on gauge invariant states is naturally implemented. 

Discretization of the Dirac operator is analogous to the now standard 
construction of the supersymmetric Yang-Mills theories on a lattice ||. We 
do not address here important, and specific for D=10, questions of Euclidean 
formulation for the fermionic degrees of freedom already discussed in [H and 



Weyl projection on the lattice |10|- Due to the Ni — 1 periodicity all hop- 
ping terms along the space directions collapse into the diagonal blocks of the 
fermionic matrix. Then one faces the problem of evaluating fermionic deter- 
minant or pfaffian for Majorana constraint. For one dimensional system the 
fermionic matrix is effectively block three diagonal. This allows for impor- 
tant numerical simplifications. For example, we have developed an algorithm 
which reduces the computational effort of the exact evaluation of the pfaf- 
fian of the antisymmetric fermionic matrix from 0(V 3 ) to 0(V), V being the 
volume of the system. Even with this improvement, however, lattice simula- 
tions with dynamical fermions are much more time consuming than the pure 
gauge computations. 
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On the other hand experience in lattice QCD shows that the effect of dy- 
namical fermions is mostly accounted for by a shift of the coupling constant 



while neglecting the functional dependence of the fermionic determinant fllT 
Therefore, our first simplifications is to study (Q) in the quenched approx- 
imation. Indeed, as is seen below some essential features (e.g. the general 
phase structure) of the model are preserved. 

Next simplification concerns the large N limit. Again, results of lattice 
simulations for QCD strongly suggest that 't Hooft choice of the coupling 
constant, A = g 2 N, takes into account main large N effects. In fact it 
was found that even results for SU(2) gauge group are not far from those 



with higher N |12|]. We therefore propose to study systematically the N 
dependence of some features of the model (|2J), beginning with N = 2. 

Finally, present simulations are done for D — 4. Although the ultimate 
goal is D = 10, we expect that the quenched approximation, which we study 
numerically here, behaves smoothly with D []. Although this will be differ- 
ent in the full unquenched simulation, we have decided to study first the 
properties of the simpler system with D = 4. 

Needless to say, one should gradually remove the above approximations 
in the forthcoming computations. This is especially important for studying 
the low temperature phase where the supersymmetry restoration may be 
essential. 

2. Results. One of the most exciting feature of the new theory is the expla- 
nation of the Bekenstein-Hawking entropy puzzle in terms of the microscopic 
degrees of freedom of the elementary strings/branes [I4~|. In particular the 



theory predicts existence of the phase transition at which a black hole dis- 
solves into its elementary constituents [fl5|, |^|. Confronting this property with 
the predictions of the matrix model would provide an important test of the 
BFSS conjecture. Moreover, lattice study of QCD at finite temperature show 
that the very fact of the existence of the phase transition is not sensitive to 
quenching. Therefore, as a first application of our construction, we have 
chosen to study the phase structure of the system (0). Obviously, the one 
dimensional system with local interactions cannot have any phase transition 
for finite N, but just a crossover between two types of the behaviour. How- 
ever for infinite N a sharp transition may occur [[UJ. Thus the number of 



colours plays a role similar to a volume in statistical systems. Consequently, 



1 This is also confirmed by recent analytical results p 
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N t 


Plow 


Pup 


2 


1.25 


1.5 


3 


3.5 


5.0 


4 


8.0 


16.0 


5 


15.0 


40.0 


fit: 


: & = 




X 2 /NDF 


a 


7 


.55/2 


.17 ± .05 


3.02 ± .33 



Table 1: Estimated location of the transition region /3 C e (Plow, Pup) for 
different lattice sizes iV t and results of the power fit. 

we expect some signatures of the phase change for finite and even small N. 
Subsequent simulations for larger N would provide more information for the 
quantitative (e.g. finite size scaling) analysis. 

As an order parameter we choose the distribution of the Polyakov line 

p = ^ UdW) ■ ( 4 ) 

Similarly to lattice QCD, symmetric concentration of the eigenvalues around 
indicates a low temperature phase (which would have the interpretation 
of a black hole phase in the full model) where (P) ~ 0, while clustering 
around ±1 (for SU(2)) is characteristic of the high temperature (elementary 
excitations) phase. 

A sample of results for different N t (= No) and f3 is shown in Fig. 1. 
Indeed, for each Nt, we see a definite change of the shape with f3. This is 
the first result: the system (0) shows unambigously the onset of the phase 
change, even in the quenched approximation and for N = 2. 

Second, the dependence of the pseudocritical temperature P c on the time 
extent Nt is consistent with the continuum limit expectations T c ~ (g 2 N) 1 ^ 
fT3| . Indeed, the temperature of a system is given by T = l/(aN t ) . Together 
with (|3|) these relations imply P c ~ Nf. The estimates for P intervals where 
the change of phases occur are presented in Table 1 for several lattice sizes 
N t . Results of the power law fit are also quoted. A good quality of the 
fit and the agreement with the canonical exponent, 7 = 3, is encouraging. 
Simultaneously, we obtain the proportionality coefficient a which translates 
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into 

T c = (^) V VA0 1/3 = (-28 ± M)(g 2 N)^. (5) 

To summarize this point: the observed dependence of (3 C on N t agrees with the 
canonical scaling expectations for the one dimensional system, and indicates 
the finite value of the transition temperature in the continuum. Moreover, 
the coefficient in the continuum relation (JD has been determined for the 
first time. Only proportionality of the two scales has been considered until 
now |n| 0, [L3||. Since both the pseudocritical temperature and a depend 
in general on N, it is important to repeat similar analysis for higher gauge 
groups. 

Next we study the temperature dependence of the total size of the system 
R 2 = 9 2 Ea(A?) 2 0. We define for SU{2) 

(R 2 ) = (4 - ((Tr(f/ S )) 2 )) /a 2 , (6) 

where U s is any space link. Due to the periodicity (N s = 1) in space (H) is 
gauge invariant. The space links U s are the remnants of the torelon observ- 
ables well known in lattice QCD fl7| . 



One dimensional Yang-Mills coupling g provides a single scale for all 
continuum observables similarly to Aqcd in four dimensions. In the following 
all dimensional quantities quoted in units of g 2 ^ 3 are denoted by a tilde. 

Even though the quantum mechanical system (||) is much simpler than the 
full .D-dimensional field theory, extracting the continuum limit of the lattice 
formulation (|2|) may be a nontrivial task. For example, the above limit con- 
tains the complete information about both the perturbative weak coupling 
and nonperturbative strong coupling regimes in the continuum. Technically, 
relation @ implies that a reasonably small lattice constant, a requires simu- 
lation with a very large coupling (3. In addition, the one dimensional systems 
are harder to thermalize. All this poses an interesting challenge in construct- 
ing new algorithms suitable for this problem. Some of such alghorithms are 
under development and will be discussed elsewhere. Here we use mostly the 
standard local Metropolis update. To overcome the critical slowing down 
we simply increase the number of thermalization and decorrelation sweeps 
with (3, until results become independent of the starting configuration. This 
turned out to be in accord with the dynamical exponent z — 2. For example 
when running at a = 1.0 we used 5000 thermalization and 50 decorrelation 
sweeps, while for a = 0.1 about 10 6 thermalization and 5000 decorrelation 
sweeps were required. One of the new algorithms mentioned above is the 
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SU(2) heat bath designed for an update of the space-space plaquettes in 
(FJ) which contain twice the same link. Current version is effective only for 
(3 < 64. Results obtained with the new heat bath and the standard Metropo- 
lis agree within statistical errors. To check independently the performance of 
the Metropolis algorithm for higher f3 we have also monitored the correlation 
length in the torelon channel at zero {i.e. low) temperature. It reveals the 
expected canonical scaling with a. 

Fig. 2 shows the dependence of R 2 on a, for several values of the temper- 
ature T. MC results depend smoothly on a, at fixed T, which confirms the 
existence of the continuum limit ([]). The a dependence is clearly different 
in low and high temperature regions. For .1 < T < .3, (R 2 ^ is practically 

independent of a and points for different (but small) T collapse on the same 
line. For higher T quadratic minimum at a = develops and shrinks with 
the further increase of the temperature. For T > 1.5 simulations for smaller 
a are required in order to see this structure and determine the continuum 
limit. We have also extracted (R 2 ) from another lattice observable \Tr(U s )\ 
with practically the same results. 

Fig. 3 shows the size of a system extrapolated to a = as a function 
of the temperature. Both quadratic and quartic fits of a dependence were 
used to perform the extrapolation [|18f| . We have also checked the stability of 
quadratic fits with respect to removing one or two data points with smallest 
a (highest (3). Results of the extrapolation were stable with respect to all 
these variations. Small systematic shifts are included in the errors displayed 
in Fig. 3. The location of the transition region is in rough agreement with the 
estimate (||) of the pseudocritical temperature T c = 0.35 ±0.04 0. Again, it is 
evident that the system is indeed different in the two regimes. Moreover, our 
results agree qualitatively with the analytical prediction obtained by solving 
a gap equation in the infinite N limit [O]. The latter gives a temperature 



independent constant at low temperatures and the classical T 1//2 growth for 
high temperatures. We have also found a reasonable agreement with a sim- 
ple mean field model for SU(2) with the gauge projection^. As expected 
the model does not have a phase transition, but shows a smooth crossover 
located as in Fig.3. The constant value for (R 2 ) is satisfactorily reproduced 
in the low temperature, vacuum driven region. At higher temperatures the 



2 The pseudocritical temperatures determined from different observables can be differ- 
ent. 

3 To be discussed in detail elsewhere. 
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model predicts intermediate linear, albeit weaker than MC, behaviour which 
asymptotically turns over into T 1//2 as in the infinite N case. 
3. Conclusions. We have constructed the matrix model of M-theory on a 
lattice in D=2,4,Q and 10 dimensions. The resulting system corresponds 
to the supersymmetric formulation of Yang-Mills theory on the asymmetric 
-D- dimensional lattice with all D — 1 space extensions N s — 1. The new con- 
struction was tested in the quenched approximation for D=4. In particular, 
we have found the onset of a black hole to strings transition even for the 
SU(2) gauge group. The pseudocritical temperature was determined. The 
size of the system was also measured at different temperatures and lattice 
cut-offs. It shows the expected canonical scaling. After extrapolation to 
the continuum limit it confirms the existence of the two phases and agrees 
qualitatively with the mean field calculations. 

A host of new applications can follow. On the technical side, new algo- 
rithms are required to reduce the critical slowing down at very large values 
of the lattice coupling. Such studies have already begun. Including dy- 
namical fermions is facilitated by the linear nature of the system and may 
lead to more efficient fermionic algorithms. Certainly the issue of dynami- 
cal fermions is very important especially in the low temperature phase since 
one expects that supersymmetry should be broken only in a minimal fashion 
there. With dynamical fermions in D=10 one may have to use the recently 



proposed chiral formulation On the other hand for the reduced system 
the task may be simpler than e.g. for QCD. It would also be very interesting 
to apply analytical methods developed in |19|, |20|. Incidentally, a merit of 



the present approach is the possibility to draw from the expertise, techniques 
and algorithms developed in the lattice community. 

A systematic study of the model for higher N would allow finite size 
analysis and determine more detailed characteristics of the transition. In 
particular it would be interesting to check if the "soft" dependence on N 
observed for 3D and 4D SU (N) lattice YM fI2|| , persists in the SYM quantum 
mechanical model. 

Finally, one of the ultimate physical goals would be to study the ther- 
modynamics of the black-hole phase in the full D=10 model and verify exis- 
tence of the rich phase structure predicted by the string/M theory ||. This 
would also provide a possible nontrivial quantitative test of (a version of) 
the AdS / CFT correspondence at strong coupling not protected by any non- 
renormalization theorems P, |2T|. 

Last but not least, many other problems inspired by the BFSS conjecture 
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can be quantitatively studied. 
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Figure 1: Distribution of the Polyakov line (f|), — 1 < P < 1, for different (3 
and N t . Note different /3 range for different N t . 
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Figure 2: Dependence of the total size of the system on a, for T = .1 — 
.3, .6, .9, 1.2 and 1.5 (upwards) in units of g 2 ^ 3 . Quartic fits are represented 
by the solid lines. 
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Figure 3: Size of the system (§) extrapolated to the continuum, as a function 
of the temperature. 
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